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Abstract 

We prove that the variance of the current across a characteristic is of 
order t 2 ^ 3 in a stationary asymmetric simple exclusion process, and that 
the diffusivity has order t 1 ^ 3 . The proof proceeds via couplings to show the 
corresponding moment bounds for a second class particle. 
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1 Introduction 

The asymmetric simple exclusion process (ASEP) is a Markov process that de- 
scribes the motion of particles on the one-dimensional integer lattice Z, subject to 
the exclusion interaction that allows at most one particle at each site. Particles 
in the process jump one step to the right with rate p and one step to the left 
with rate q = 1 — p, and we assume < q < p < 1. Particles attempt jumps 
independently of each other, but any attempt to jump onto an already occupied 
site is suppressed. In Section [2] below we give a rigorous construction of ASEP in 
terms of Poisson clocks that govern the jump attempts. This process is among 
the interacting particle systems introduced in Spitzer's seminal paper [25]. We 
refer the reader to Liggett's monographs [THE] for coverage of most of the work 
on ASEP up to the late 1990's. 

In 1994 Ferrari and Fontes [TO] proved a central limit theorem for the net 
particle current seen by an observer moving at a fixed speed v. This quantity 
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that we denote by J^"'(t) is the number of particles that pass the observer from 
left to right minus the number that pass from right to left during time interval 
(0, t]. The particle process is assumed to be stationary with Bernoulli-distributed 
occupation variables at some average density g 6 (0, 1). The result is the weak 
limit of the diffusively rescaled and centered current: 

,. jW(i)-E[jW(t)l 
, hm 7i72 = 

The limit \v is a centered Gaussian random variable with variance 

<? 2 = g(l-g)\(p-q)(l-2g)-v\. 

The interesting phenomenon is the vanishing of the variance at the character- 
istic speed v = V e = (p — q)(l — 2g). As we explain below, V s is the speed at 
which perturbations travel in the system, both at the microscopic particle level 
(in the expected sense as given in (|2.7[) below) and at the macroscopic p.d.e. level. 

Physical reasoning [5^ implied that the correct order of the fluctuations of 
the current across the characteristic should be i 1 / 3 . This would explain the 
degenerate limit under i 1 / 2 normalization. These "i 1 / 3 fluctuations" remained 
elusive throughout the 1990's. 

The seminal papers of Baik, Deift and Johansson [5] and Johansson [T3] gave 
the first rigorous proofs of such fluctuations. The correct order was verified to be 
i 1 / 3 , and the limiting fluctuations were found to obey Tracy- Widom distributions 
from random matrix theory. The first paper dealt with the last-passage version 
of Hammersley's process, and the second with the last-passage version of the 
totally asymmetric simple exclusion process (TASEP). Total asymmetry means 
here that particles are allowed to jump only in one direction at a constant rate, 
so this is the case p = 1, q = 0. 

These papers did not study stationary particle processes, but instead processes 
started from special jam-type deterministic initial conditions. For TASEP this 
means that initially all sites to the left of the origin are occupied and all sites to 
the right empty. With such initial conditions the processes could be represented 
by last-passage percolation models, a point that had been exploited already in 
the past (among the seminal ones were [T], [22] and [23]). The actual analysis was 
then performed entirely on combinatorial descriptions of the last-passage model. 
Later a last-passage representation was also found for a stationary TASEP [18], 
and then the Tracy- Widom limit proved for the current across the characteristic 
in that setting [TJ] . 

The early proofs of fluctuations relied on a counting argument that utilizes 
the Robinson-Schensted-Knuth correspondence for Young tableaux, and Gessel's 
formula that converts certain Schur function sums into Toeplitz determinants. 
Later this step has been replaced by a more direct connection between the last- 
passage model and a determinantal point process. The fluctuation limits are 
then derived by analyzing the asymptotics of the determinant in the appropriate 
scaling regime. Consequently, while a genuine breakthrough has been achieved, 
the delicate steps of the proof restrict the reach of the results in several ways. 
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In particular, the particles of the systems are permitted to move in only one 
direction and admit only the simplest type of jumps. 

In the present paper we give the first proof of the accurate order of the fluc- 
tuations in systems that are only partially asymmetric. Namely, we show in the 
original setting of Ferrari-Fontes [TU] that the variance of the current across the 
characteristic in (p,q) ASEP is of order i 2 / 3 . Our arguments are entirely prob- 
abilistic and utilize couplings of several processes and bounds on second class 
particles. Informally speaking, second class particles are perturbations in the 
system that do not disturb the motion of the regular particles but are influenced 
by the ambient system. Precise definition must wait for the construction of the 
coupled processes in Section [21 Presently it appears that there is no way to apply 
the combinatorial-analytic approach pioneered in [2] to ASEP because there is 
no last-passage model where the analysis could begin. 

We take a key insight from recent work of Cator and Groeneboom [7j and 
from our joint work with Cator in [4|: this is the idea of coupling processes 
whose densities differ by 0{t~ 1 / z ) in order to bound the motion of a second class 
particle, whose fluctuations in turn are linked to the variance of the current. The 
couplings we utilize go back to Ferrari, Kipnis and Saada [11] who introduced 
them to study the microscopic locations of shocks in the particle system. 

Fluctuation results for asymmetric exclusion processes have also been stated 
in terms of a quantity called the diffusivity D(t). One way to view the link 
between current variance and diffusivity involves the second class particle: the 
variance of the current is the expected absolute deviation of the second class 
particle, while tD(t) is the variance of the second class particle. For ASEP we 
also obtain the correct order t 1/>3 for the diffusivity. In Section [2] we state the 
main result Theorem 12.21 which is a moment bound for the second class particle. 
Bounds for current variance and diffusivity appear as Corollaries 12.31 and 12.41 

There is also work on the diffusivity with resolvent methods. The resolvent ap- 
proach can handle more general jump kernels than the nearest-neighbor type (the 
modifier "simple" in ASEP refers to the restriction to nearest-neighbor jumps). 
But so far this approach has not yielded optimal bounds. Results in one and two 
dimensions have appeared in [T5] and [27] . 

The most recent work utilizing the resolvent approach is by Quastel and Valko 
[20l [T9] . In [20] they show that, for any two finite-range exclusion processes 
with nonzero-mean jump distributions, the ratio of Laplace transforms of tD(t) 
is uniformly bounded. In 19 this comparison theorem is paired up with our 
Corollarv l2.4l to get the correct order of the Laplace transform oitD(t) for all these 
exclusion processes, and with an additional argument also the correct pointwise 
upper bound D(t) < Ci 1 / 3 . 

A few more words about the broader context. We see this paper as an open- 
ing for a treatment of several other models, such as zero-range and bricklayer 
processes. For these systems diffusive current fluctuations off the characteristic 
were established by Balazs in [3J. 

The t 1 / 3 fluctuations with Tracy- Widom limits are universal for some class of 
asymmetric systems whose precise characterization is not clear at the moment. 
There is also another universality class among asymmetric systems in one di- 
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mension, one where fluctuations occur on the scale i 1 / 4 and limits are Gaussian 
processes related to fractional Brownian motion with Hurst parameter H = 1/4. 
Such results appear in the papers [5], [53], and [S]. 

Notation. We adhere to the usual notation for particle systems, except that 
we underline notions that pertain to the entire lattice: so rj(t) = {??i(£)}iez 
denotes the occupation variables of an exclusion process, and /i = /i® z is the 

i.i.d. product measure with marginal fi. In general X denotes a centered random 
variable: X = X — EX. C and Cj for i = 1, 2, 3, . . . denote positive constants 
that can depend on g and p and can change from line to line. 

Acknowledgement. We thank an anonymous referee for a thorough reading of 
the manuscript and insightful comments. 

2 The exclusion process and the results 

Construction of the process and second class particles 

The asymmetric simple exclusion process (ASEP) is a Markov process on the 
state space f2 = {0, 1} Z . Given a state w = {u>i}i^z S the following jumps can 
happen independently at different sites: 

(2.1) (uji, u i+l ) — > (cji - 1, Lu i+ i + 1) with rate pui(l - oJi+i), 

(2.2) (u>i, uJi+i) — ► (wi + 1, uJi+i - 1) with rate q(l - W;)a;j+i. 

We assume < q = 1 — p < p < 1. The special case p = 1 is called TASEP, or 
the totally asymmetric simple exclusion process. 

We interpret the process as representing unlabeled particles that execute inde- 
pendent nearest-neighbor random walks on Z, subject to the exclusion interaction 
that suppresses attempts to jump to an already occupied site, p is the rate of 
a particle to jump to the right and q is the rate to jump to the left. The value 
u>i(t) = means that site i is vacant at time t, and u>i{t) — 1 that site i is occupied 
at time t. The state of the entire process at time t is then uj(t). 

A rigorous construction of this process is done by giving each site i two Poisson 
processes on the time line [0,oo): a rate p process A^_>i+i and a rate q process 
Ni->i-\. The processes {JVj_j + i, Ni-^-i :ieZ} are mutually independent, and 
also independent of the initial configuration w(0). The rule of evolution is that 
when Ni^i + \ jumps, a particle is moved from i to i + 1 if i is occupied and i + 1 
is vacant. And similarly with iVj_,i_i. Thus the rates (|2.ip ~ (|2.2p are realized. 

Let n g denote the measure /i e {l} = g = 1 — ^^{0} on the set {0, 1}, and let 

= ^ e ^ e i-i-d- Bernoulli product measure with mean density g on fi. It is 
known that the measures {/i : < g < 1} are the extreme points of the convex 
and weakly compact set of invariant distributions for the process that are also 
invariant under spatial translations. 

It is convenient to embed the exclusion process in a height process that rep- 
resents a wall of adjacent columns of bricks. On top of each interval [i, i + 1] 



4 



sits a column of bricks with height hi G Z. The entire height configuration is 
h = {hi}i e z, restricted to satisfy 

(2.3) < hi-i — hi < 1 for each i 

so that the wall slopes downward to the right. Let the Poisson processes govern 
the evolution of the heights: when JVj—vi+i jumps add a brick on top of the column 
on [i, i + 1], and when Ni+i—>i jumps remove a brick from the column on [i, i+ 1], 
But suppress every step that leads to a violation of (|2.3p . 

Given an initial particle configuration lo(0) , define an initial height configura- 
tion by 

fE°=iW°) fori<0, 
h l (0)=lo fori = 0, 

[-Ei=i^(0) fori>0. 
Let the heights evolve, and define 

uji(t) = hi-\(t) - hi{t). 

Then this process tj(t) is exactly the ASEP constructed earlier, and the height 
increment hi(t) — hi(0) is the net particle current across the bond (i, i + 1). 

The Poisson construction reveals its power when it is used to run simulta- 
neously several processes started from different initial states. This is called the 
basic coupling. The first observation is that this coupling preserves monotonicity 
among both particle and height configurations. Ordering is defined sitewise: for 
particle configurations i] < lo means that r\n < uji for each i 6 Z, and similarly for 
height configurations g < ft, if gi < hi for each i 6 Z. The basic coupling has this 
property, called attractivity: 

77(0) < w(0) tjO) < u(t) and g(0) < h(0) g(t) < h(t) 

for all t > 0. 

We use the following terminology: if we have two coupled exclusion processes 
such that rj(t) < io(t), then the cu — 77 second class particles are the particles 
that occupy sites i at which u>i(t) — r)i(t) = 1. The joint process (j](-),ui(-)) can 
be constructed from a two-class process: (i) The first class particles 77 obey the 
ASEP dynamics as described earlier, (ii) The second class particles di — uji — r\i 
also obey the Poisson clocks when they can, but they are not allowed to jump on 
sites occupied by first class particles, and when a first class particle jumps on a 
site occupied by a second class particle, they swap sites. 

Let 5$ £ fi denote a configuration that has only a single particle at site i. 
If r) E Q is such that 770 = 0, we can legitimately define r/ + = i] + S Q . In this 
situation we say that there is a single second class particle between ?y + and T] at 
site 0. Since the basic coupling conserves the single second class particle, there 
is always a site Q(t) such that 

(2-4) V + (t)=r[(t)+S Q{t y 



5 



Q(t) is the position of the second class particle at time t, which performs a nearest 
neighbor walk, influenced by the ambient process rj(-). 

It is convenient to also have the notion of a second class cmfiparticle at po- 
sition Q a (t) in a process to(t). This means that Q a (t) is the location of the 
single discrepancy between two processes and uT (t) that are started so that 
w~(0) = lu(0) — where i = Q a (0). A moment's reflection reveals that in fact 
in the basic coupling of ASEP a second class particle and an antiparticle are the 
same thing. But in the proofs we will couple more than two processes and this 
extra flexibility will be convenient. 

Current fluctuations and diffusivity 

Let [x] denote the first integer from x towards the origin, in other words [a;] = [^J 
(floor) when x > and [x] = \x~\ (ceiling) when x < 0. For a speed value V 6 R 
define 

(2.5) J^(t) = h [vt] (t), 

the height of the column over interval [[Vt], [Vt] + 1] at time t. Due to the 
normalization ho(0) = 0, j( v \i) is the total net particle current seen by an 
observer moving at speed V during time interval [0, t]. Or more concretely, 
j( v Ht) = j| V) (t) - J (y \t) where J { +\t) is the number of particles that be- 
gan in (— oo,0] at time but lie in [[Vt] + l,oo) at time t, and J^\t) is the 
number of particles that began in [1, oo) at time but lie in (— oo, [Vt] ] at time 
t. One can compute the density g equilibrium expectation 

Ej( v Ht)=t(p-q)g(l-g)-g[Vt] 

by writing a martingale for ho(t) and then adding in h[ V t] (*) — h (t) which counts 
particles between sites and [Vt]. 

Our results are based on an interplay between currents and second class par- 
ticles. One key fact is the next connection. 

Proposition 2.1. Let lu(-) be an ASEP started from its stationary Bernoulli 
distribution /i^. Condition the origin to be empty and let Q(-) be a second class 

particle that starts at Q(0) — 0. Alternately condition the origin to be occupied 
and start a second class antiparticle Q a (-) at the origin. Either situation can 
be used to compute the variance of the current of the stationary process for any 
V € M: 

Var(jW(t)) = g(l g)V{ \ [Vt] Q(t)\ \ c (0) = 0) 

= e(i-e)E(|[Vt]-Q„(t)||wo(o) = i) 

We shall also have occasion to use the following identity (true under the 
assumptions of the proposition above): 

(2.7) E(Q(i) | c^o(0) = 0) - E(Q„(i) | c^ (0) = l) = t(p - q)(l - 2g). 
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Equations (|2.6p and (|2.7p have been derived and utilized earlier for TASEP, 
by Ferrari and Fontes [10] and (|2.7|) by Prahofer and Spohn [18]. In article [6] 
identities (|2.6p and (|2.7p are proved not only for ASEP but also for zero-range 
and bricklayer processes. 

The interesting current fluctuations occur at the characteristic speed V e — 
(p — q)(l — 2g). From (|2.7p we see that this is the average speed of the second 
class particle. To explain the significance of V e from another perspective, we 
digress briefly to discuss a large scale, deterministic view of exclusion dynamics. 

On large space and time scales the macroscopic particle density g(t 7 x) of 
ASEP obeys the conservation law 

(2.8) g t + f(g) x = 

with flux f(g) = (p — q)g(l — g). This notion is made precise through a law of 
large numbers called a hydrodynamic limit. One considers a sequence of processes 
a/ n ) (t) such that, for each n, the initial occupation variables are independent with 
means Ea4 n } (0) = eo(i/n) for a given continuous function g . Then we have the 
limit 

b 

(2.9) lim — cul n \nt) — / g(t,x) dx in probability 

na<2<n& „ 

for all intervals [a, b}. The limiting density g(t, x) is the entropy solution of (|2.8|) 
with initial condition g(0, x) = go(x). Away from shocks the characteristics of the 
partial differential equation (|2.8|) are solutions of the ordinary differential equation 
x = f'(g(t,x)). At constant density g the characteristic speed is f'{g) = V s - 

Limit (|2.9p has been proved over time by a number of authors under various 
hypotheses. Article [21] was the first to achieve a satisfactory level of generality 
for a broad class of asymmetric systems, also in higher dimensions. For more in- 
formation we refer to [14] for hydrodynamic limits of interacting particle systems 
and to [9 for the theory of scalar conservation laws. 

Let us return to the stationary ASEP with Bernoulli-^ occupation variables 
at each fixed time. A basic object for understanding space-time correlations is 
the two point function S(i,t) = E[wi(i)u;o(0)] — g 2 . It can be written as the 
transition probability of the second class particle: 



(2.10) 



S(i, t) = g(E[wi(t) | wo(0) = 1] - E[wi(t)] ) 

= g(l - e)(E[«i(t) | wo(0) = 1] - E[ui(t) | Wb(0) - 0] ] 
= g(l - e)E[ Wi + (t) -w,(t) |w+(0) = 1, wb(0) = 0] 

= e(i-e)P(Q(t) = i|wb(o) = o). 



On the third line above we coupled w(-) with a process w + (-) whose initial con- 
figurations satisfy w + (0) = w(0) + S Q . In other words, the second class particle 
resides initially at the origin. Note that use of the term "transition probability" 
is not meant to suggest that Q(t) is a Markov process. 
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From (|2.10p and (|2"77|) follow the properties 

J2S{i,t) = g(l- g) and ^iS(i,t) = g(l - q)VH. 

The diffusivity is by definition a normalized second moment of the two-point 
function: 

(2.ii) ^) = ^)E(-^(^)- 

Consequently the diffusivity can also be expressed in terms of the variance of the 
second class particle: 

D(t) =t- 1 Var{Q{t)\u J o{0) = 0). 

If the second class particle behaved like a random walk (diffusively) then its 
variance would be of order t, and thereby D(t) of constant order. The interesting 
point is that the ASEP second class particle is superdiffusive with variance of 
order i 4 / 3 . 

We can now state the main theorem, a moment bound for the second class 
particle in the setting where all the other occupation variables start off Bernoulli- 
distributed. 

Theorem 2.2. For ASEP w(-) with rates < q = 1 — p < p < 1, started in \i 
distribution with any g G (0, 1), there exist constants < to, C < oo such that for 
real 1 < m < 3 and t > to 



(2.12) C" 1 < E 



Q(t) - VH 

t 2/3 



w (0) = \ < c 



Let us simplify notation to J e [t) — J^ ye \t) for the current across the char- 
acteristic. Via (|2.6p . taking m = 1 gives the variance of this current. 



Corollary 2.3. With assumptions as in Theorem \2.SX fort > to the variance of 
the current across the characteristic satisfies 

C -l t 2/3 < Var(J e (t)) < Ct 2/3 . 
Taking m = 2 identifies the order of the diffusivity. 
Corollary 2.4. With assumptions as in Theorem Y&.'A for t > to 

c-h 1 ' 3 < D(t) < ct x '\ 

The upper bound of (|2.12p is not valid for alH > 0. For small t the variance 
of Q(t) is of order t because the second class particle is likely to have experienced 
at most one jump. 

The rest of the paper is devoted to the proof of Theorem 12 . 21 first the upper 
bound and then the lower bound. 
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3 Upper bound 



Proof of the upper bound utilizes couplings of several processes and a tagged 
second class particle introduced by Ferrari, Kipnis and Saada [11] . We refer the 
reader to the exposition of these ideas in Section III. 2 in Liggett's second mono- 
graph [17) . Before deriving the upper bound we introduce some preliminaries on 
the couplings. 

Couplings 

We start by describing an initial distribution on three classes of particles and 
with labels attached to particles of the two lower classes. 

Fix densities < A < q < 1. Define the measure \x on {0, 1} x {0, 1} by 

(3.1) fi(0,Q) = l- g, fi(0,l) = g-X, /x(l,0) = 0, and /j(1, 1) = A. 

The first marginal of fi is n\ and the second is /i e . Let fi = be the i.i.d. 
product measure on [{0, 1} x {0, 1}] Z with marginals fi. (Our notation for /j, 
differs from that in [T7] where the second variable of fi is the process of second 
class particles.) We condition this measure on having an w — r\ second class 
particle at the origin, and define the measure p, by 

fi{ d?7, dui) = /_t( dr?, du | r/ a = 0, ui = 1). 

Then /2-a.s. we have the following conditions for r\ and lu: 

• -q <u, 

• there is a second class particle between rj and w at the origin, and there are 
infinitely many of them on both sides of the origin. 

Label the lu — r\ second class particles with integers in an increasing fashion from 
left to right, giving label to the one initially at the origin. Xk(Q) denotes the 
position of second class particle with label fc, so that 

(3.2) • • • < X_i(0) < AT o (0) = < Xi(0) < • • • 

Given (77, lu) chosen from distribution /x, define a third configuration £ as 
follows. Let = uji whenever r\i = lui. At all other sites i we have an lu — r/ 
second class particle: = 77, < u>i = 1. If this second class particle has label n, 
then let 



i]i = 0, with probability 



( 3 - 3 ) (1 ~ Cx„(o) 



LUi = \, with probability 



1 + (p/q) n ' 

ip/q) n 
1 + (p/g)« ' 



where these choices are independent from each other and everything else. For 
the case p = 1 the probabilities are defined as the limits of the above formulas as 
V =\-q/\. 
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Let Jl* denote the resulting joint distribution of (77, £, uf). It is clear that /l*-a.s. 
we have 77 < £ < uj. 

Let these configurations evolve from the initial distribution Jl* with common 
Poisson clocks. This basic coupling preserves the ordering n{t) < ((t) < uj_(t) a.s. 
for all times t. The effect of the coupling is that the C ~ V particles have priority 
over the ui — £ particles. 

Consider the evolution of the the uj — rj second class particles in this coupling. 
They start off labeled with integers as described above (|3.2j) . We let each particle 
keep its integer label and let Xk (t) denote the position of particle labeled k at time 
t. Since only nearest-neighbor jumps are permitted, the ordering Xk(t) < Xk+\{t) 
is preserved. Additionally, particle k carries a £-mark Cx k (t)(f) G {Oj 1}> initially 
assigned by (|3.3p . Cx k (t)(t) = indicates that in the further subdivision into 
second and third class particles, at time t position Xkit) is occupied by a third 
class particle. These £-marks are not retained but exchanged. Whenever Xk 
and Xk+i are neighbors and carry different £-marks (1,0) or (0,1), then (1,0) 
becomes (0, 1) with rate p and (0, 1) becomes (1,0) with rate q. This is just a 
restatement of the effect of the basic coupling when a £ — r\ particle is next to an 
u) — C particle. 

The oo — r) second class particle that starts at the origin is of special importance 
to us, so we set X{t) = X (t). 

Proposition 3.1. Couple the processes as described above with initial distribu- 
tion Jl* . Then for any fixed time t 6 [0, 00) the marks {Cx k (t)(t) : k £ Z} are 
independent of the entire process {(rj(s) , uj_(s)) : s > 0}. Furthermore, the dis- 
tribution of marks is constant in time: given the evolution (t)(-), for all 
< t < 00 and disjoint sets Iq and I\ of integer labels, 

p {Cx k (t)(t) = for k e Io and Cx n (t)(t) = 1 for n & h \ ??(•), w(-)} 

= tt i . n {p/q)n 

11 ! + {p/q)k 11 ! + {p/q) n 



This statement is hinted at in a remark after Proposition III. 2. 13 in [T7] , 
We summarize the argument for the reader's convenience. Observe first that we 
can construct the process of (?7,a;)-particles and ^-labels with two independent 
collections of Poisson clocks, one that governs the process (rj(-) , , and another 
one that governs the exchanges of (- marks among the to — rj second class particles. 
This construction is not the same as the basic coupling described above, but it 
leads to the same process because the infinitesimal rates are the same. 

Consequently we can first construct the (j](-),oj(-)) evolution for all time, and 
then superimpose on it the evolution of the £-marks. The claim is that, given 
the entire evolution (??(•), w(-)), the mark configuration {(x k (t)(t) : fc e Z} on the 
uj — rj second class particles has distribution (|3.3p . 

The key point is reversibility. Let us abbreviate temporarily Uk (t) — Cx k (t) (t) 
for the mark of second class particle Xk at time t. If the process w(-) obeyed 
the (p, <7)-ASEP dynamics, the product distribution with marginals (|3.3|) would 
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be one of the well- known blocking measures that are reversible for u(-). (This 
would in fact be the case if A — and g = 1, for then the second class particles 
would occupy all sites and never move: -Xfc(i) = k for all < t < oo.) The result 
now follows from this observation: the evolution (??(•), ^(-)) can be thought of 
as a dynamical "environment" for the mark process that admits or prohibits 
certain exchanges at different times: X}. and Xk+i cannot exchange marks unless 
they occupy adjacent sites, and the time intervals during which this happens 
are determined by (tj(-),uj(-)). However, imposing such an environment on the 
process does not change the reversibility of the measure. 

This last point can be checked rigorously by first letting only finitely many 
marks Uk^) = {v,k(t) : —K < k < K} evolve while the remaining marks are 
frozen. Then we are talking about a finite state Markov chain. Given the evolu- 
tion (j) (■),&(■)), we can partition the time axis = to < t\ < t% < ■ ■ ■ < U / oo 
so that the generator of the chain Uk(-) does not change during (ti-x,U). Since 
detailed balance is not violated by symmetrically prohibiting certain jumps, re- 
versibility and hence invariance is preserved during each interval (tj_i,tj), and 
thereby for all time. Letting K — > oo extends the invariance to the infinite mark 
process. This argument establishes Proposition 13. II 

The purpose of the construction is to confine certain special particles to be 
introduced shortly. Notice that a.s. there are only finitely many positive labels 
n > for which the first event of (|3.3[) happens. Hence we have a rightmost 
position R(0) — sup{z : r\i — Q < LUi}. This is the initial position of the rightmost 
second class particle between £ and u. Similarly, the second event in (|3.3|) happens 
a.s. for only finitely many negative n < 0. Let L(0) = mf{i : rji < Q = uii} be the 
leftmost position where this happens, the initial position of the leftmost second 
class particle between rj and £■ Define the events 

A = {0 < i?(0)} and B = {L(0) < < R(0)}. 

These events only depend on the marks (|3.3p . and hence are independent of (rj, lu). 
A always has positive probability, while B has positive probability if p < 1. 

These positions evolve in the coupling, so at time t we let R(t) be the position 
of the rightmost second class particle between and and similarly L(t) 
is the position of the leftmost second class particle between r](t) and ((t). L(t) 
and R(t) are always among the ui — 7/ second class particles {Xk(t)}. Let n^i) 
and nn(t) be the random labels such that L(t) = X nL u\{t) and R{t) — X nR ^(t). 
These labels are functions of the marks: 

(3.4) n £ (t)=inf{fc:C Xfc(t) (t)=l} and n fl (t) = sup{fc :&*(*)(*) = 0}. 

By the above proposition riL{t) and nn(t) are independent of (T](t), w(t)) and 
their distribution does not change with time. 

To complete the construction we define two more initial configurations by 

r] + = 7] + and uf~ = uj_ — 5 
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where S denotes a configuration with a single particle at the origin, as explained 
in Section [2] The basic coupling then applies to all five processes rj(-), t] + (-), £(•), 
lo^(-), w(-) with the above initial data. The orderings 

(3.5) 2 (t) < C(i) < w(f ) , 2 (t) < uT (t) < w(t) , ^(t) < ry+ (f ) < u#) 

are preserved by the evolution. There is always one difference between rj(-) and 
rj + (-) which is thought of as a single second class particle on ry(-) in the sense (12.41) . 
Denote its position by Q(t) with Q(0) = 0. Also, cj(-) and always have a 

single second class particle between them, thought of as a second class antiparticle 
on cj(-) as described in Section^ Its position is Q a (t), and Q a (0) = 0. See Figure 
[T]for an illustration of the initial configurations in this coupling. 
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Figure 1: A realization of the initial state of the five coupled processes. The top portion 
of the figure shows the height functions and the bottom portion the corresponding 
particle configurations. The lattice Z runs at the bottom. Solid disks denote particles 
and open circles denote vacant sites. Above the particles, * marks locations of u> — £ 
particles (first case in (|3.3f) ) and * marks locations of £ — r\ particles (second case in 
(I3.3|) h Together these make up the ui — rj second class particles and the numbers above 
these symbols represent the labeling (|3.2|l . Both events A and B happen in this picture. 
The particle-hole pairs inside the ovals represent the initial locations of Q a and Q. 

The bounds will be proved by controlling the particles Q(t) and Q a (t). The 
next lemma contains one of the key points. 

Lemma 3.2. In the coupling of five processes we have these implications: 
(3.6) A C {Q a {t) < R(t) for all t > 0}. 
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and 

(3.7) B C {Q a (t) < R(t) and L(t) < Q(t) for all t > 0}. 



Proof. We prove (|3.6[) and leave the similar argument for (|3.7p to the reader. 

Under the event A, Q a (0) = < R(0) holds initially. That is, we have aw-( 
second class particle at or to the right of the second class antiparticle Q a . We 
show that the coupling preserves the inequality Q a {t) < R(t). We have two cases 
to consider that could potentially allow Q a > R to happen. 

• Suppose the ordering £(s) < w~(s) < lo(s) holds at some time s. Since the 
ordering is then preserved for all later times t> s, CQ a (t)(t) < (t)(*) < 
w Qa(t) (*)> w hi cn implies anw-( second class particle at location Q a (t). As 
R(t) is the position of the rightmost such second class particle, Q a (t) < R(t) 
holds. 

• Consequently the only possibility for producing Q a > R is to have a jump 
in a situation of this type: Q a {t) — i, R{t) = i + 1 for some site i, and 
ur(t) = 0, Ci(t) = uj.it) = 1, Ci+i(*) = 0, w^ft) = oj l+1 (t) = 1. In this 
case column i of £ increases by one with rate p, or column i of w _ decreases 
by one with rate q. Neither one of these steps interchanges Q a and R. 

This proves (|3.6p . □ 

For applications of this lemma it is crucial that the events A and B depend 
only on the initial marks (|3.3p which were independent of the initial configu- 
ration (77, lo). Consequently A and B are independent of the joint evolution 
^(')' Q(')i When we condition the processes on A or B we 

call this construction a conditional coupling. 

Proof of the upper bound 

Throughout this section the probability measure P represents the five-process 
coupling constructed in the previous section. We begin with the proof of the 
upper bound in Theorem 12.21 C and its variants Ci,C2,... denote positive 
constants that possibly depend on p and g and whose values can change from line 
to line. We first prove a lower bound on the density of the u> — 77 second class 
particles. For integers j € Z and u > let 

j+2u 

(3.8) Nj(t)= (ui(t)-r)i(t)). 

i=j+i 

Lemma 3.3. Let A, g G (0, 1) and d > be an integer. Then there are strictly 
positive finite constants 7 = 7(g), C\ = C\{g,d) and C'2 = C2(g) such that the 
following holds: if < g — A < 7, then for all integers j G Z and u > and any 
time t > 0, 

p{Nj(t) < u{g - A) + dj < Ci exp{-C 2 u((? - A) 2 }. 
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Proof. For the moment, denote by y(-) a (p, q) exclusion process such that t/j(0) = 
uji(0) at all sites i except for i = 0, and z{-) is a (p, q) exclusion process such 
that Zi(0) = r)i(0) at all sites i except for i = 0. For i = 0, we pick the pair 
(zo(0), yo(0)) in distribution /i (|3.1[) . independently of the configuration on other 
sites. Apply the basic coupling to ensure f](t) < z(t) < y(t) < Lu{t) for all t > 
(notice that this holds initially). y(t) and z{t) are marginally time-stationary 
processes, hence we can omit the notation for their time dependence in our ar- 
guments. However, the pair (z,(t), y{t)) is not in product distribution for t > 0. 
Define 

j+2u j+2u 

Y = yj and Z = ^ Zj, 

i=j + l 



so that Nj(t) >Y - Z. Then for any a > 0, 

P{iV,(» < u(g - X) +d\ < v\c- a ^ Y - z ^ > e - Q «(e- A )- d «| 

< e au( e -\)+da _ E^ e - a ^ e Q ^) 

< e Q "(e-A)+da . [ E ( e -2oF^l/2 . [ E ( e 2aZ)]l/2 

_ e aM( e -A)+<iQ . [E(e" 2Qyo )] , ' • [E(e 2az °)] u 
= exp|au(p — A) + da + u log[l + g( e - 2Q - 1)] + u log[l + A(e 2a - 1)]} 
< exp{-au(g - A) + da + 2a 2 u(g + A - g 2 - A 2 ) + uCa 3 } 

Here we used the marginal and ^ product distributions of z and y, while the 

last inequality comes from Taylor expansion w.r.t. a. The 0(a 3 ) term is uniform 
over A < g for a fixed g. Next we pick 

g — A 



4(e + A-g 2 -A 2 )' 

this optimizes the g and A-dependent terms, and finishes the proof. □ 

We turn to develop the main estimate for the upper bound. The objective 
is to bound the deviation P{Q a (t) > u + [V e t]} with a suitable expression that 
involves the moment E\Q a (t) — [V e t}\. This is reached in (|3 . 16|) below and then 
the upper bound comes from an elementary integration step. 

Along the way we compare currents in two processes that we abbreviate as 
follows: 

J e (t) = j((p-'?)( 1 -2e))^) for currell t in the w(-) process, and 
J ve ^{t) = j«P-?)( 1 -2e))( i ) f or curr ent in the process. 

Notice that both use the same speed V e = f'(g) = (p—q)(l—2g) for the observer. 
As already defined in the Introduction, this is the characteristic speed of the w(-) 
process with density g. Let A = (P{^1}) _1 < oo from the conditional coupling, 
and recall that a tilde centers a random variable. 
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Lemma 3.4. Suppose g — X < 7 with 7 from Lemma \ 3.3l Then for positive 
integers u and times t £ [0, 00), 



P{Qa(t) >4u+[VH]} 

(3.9) < AP{JH[t) - J ye -+(t) > u(g - A) - t(p - q)(g - A) 2 } 

+ Ci exp{-C 2 u(g - A) 2 }. 



Proof. By the independence observed after the proof of Lemma [5721 conditioning 
on A does not affect the probability of the event of interest. By (|3.6[) . 

P{Q a (t) >4u+ [V e t]} = P{Q a (t) >Au+ [VH] I A} 

< P{R(t) >4u+ [VH] I A} < AP{R(t) >Au + [VH]} 

< AP{X(t) >2u+ [VH]} + AP{R(t) >Au+ [VH], X{t) <2u+ [VH]}. 

Recall that Nty et i (t) of f|3 . 8|) counts the number of oj — n second class particles at 
time t in the interval + 1, . . . , [VH] + 2u}. Since the second class particles 

stay ordered and X(t) started at the origin, the event {X(t) >2u+ [VH]} implies 
that all these second class particles crossed the path s 1— * [V e s] + 1/2 by time 
t. Each such second class particle crossing increases J B {t) — J ye ' X (t) by one. 
Therefore 

P{X(t) > 2u+[VH]} < P{J e (t) - J ve ' x {t) > N [vn] (t)} 

< P{j (t) - J ve ' x {t) > u(g - A) + 3} + P{N [vet] (t) < u(g - A) + 3}. 

Combine the previous displays to get 

(3.10) P{Q a (t) > 4u + [VH]} < AP{ J e (t) - J ye ' x (t) > u(g - A) + 3} 

(3.11) +AP{N [vet] (t) <u(q-X)+3,\ 

(3.12) + AP{R(t) > 4u + [VH], X{t) <2u+ [VH]}. 

To line (|3 . 1 1[) apply Lemma \3. 31 with d = 3. The event in (|3 . 1 2[) can happen in 
two ways: 

• The label nji(t) of R(t) is larger than u(g — A). The probability of this is, 
according to (|3.3|) and its time-invariance, bounded by 

00 00 

n=r ^-A)i 1 + ^ n= r ^-A)i p - q 

• There are fewer than u(g — A) u> — ( second class particles in the discrete 
interval {[VH] +2u,..., [VH] + 4u - 1}. The probability of this event is 
bounded in Lemma [ 
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Hence line (|3.12p is bounded by 

C l exp{-uC 2 (g - A) 2 } + Ap{p - q)- 1 exp{-uC 3 (g - A)}. 

By modifying the constants, we see that lines (|3.11[) and (|3.12[) are bounded by 
the last exponential term in (|3.9[) . 

If u and rj start from their respective [i and fj,^ equilibria, then we would 
have 

E(J«(t)) - E(J Ve ' X {t)) = Kf(Q) ~ /(A) - f(g)(g A)) 
(3.13) =-\tf"{Q){Q-\) 2 

= t{p - q)(g - A) 2 , 

where we utilized the precise form f{g) = (p — q)g(l — g) of the flux and ignored 
the error coming from the integer part of V e t. This last error is at most one. 
Our processes are also perturbed initially at the origin by the conditioning in /x, 
which gives an error not larger than 2 in the above quantity. The term +3 inside 
the probability on line (|3. 10[) makes up for these errors. □ 

Remark. We spelled out calculation (|3.13p because the nonvanishing of /" is 
a crucial factor that produces the order t 2 / 3 for the fluctuations of the second 
class particle. In the end the variable u is of order i 2 / 3 and g — A of order i^ 1 / 3 . 

Lemma 3.5. There is a constant C > such that for any time t > 0, 

E(\X{t) - L{t)\) <C{g- X)- 1 and F,(\R(t) - X(t)\) < C(g - A) -1 . 

Proof. By Proposition III. 2. 10 of Liggett |17j . the joint distribution p,S(t) of the 
pair (r)(t), as seen from X(t), can also be obtained by conditioning the 

distribution pS(t) of the pair (rj(t), at time t on having a second class 

particle at the origin. (The result is denoted by fiS(t) in [E].) The measure 
pS(t) is translation-invariant, and gives probability g — A for having a second 
class particle at the origin. If Xk(t) denotes the position of the k th second class 
particle at time t, then by Theorem B47 in [17] . we have 

E(X k+1 (t) - X k {t)) = Et^(X k+1 - X k ) = E^(X k+1 - X k ) 

= m s{t \X k+l - X k | luo - 770 = 1) 

= (f? - A)" 1 ^*) ((Xfc+i - X k ) ■ 1{^ = 1}) 
1 

g — A 

Recall that R(t) = X nii{t) (t) with n R (t) defined by (03) and X(t) = X (t). 
Use the independence and the time-invariance of the distribution of marks from 
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Proposition 13.11 

-l 

E(| J R(t)-X(t)|)=E( (Xk+i{t) - X k (t)) ■ l{n R (t) < k} 

k— — OG 

oo 

+ E( J2i X k+i(t) - X k (t)) ■ l{k < n R (t)} 

fc=0 

-l 

= E ( X k+i(t) - Xk(t)) ■ P{n R {t) < k} 



k— — oo 



J2 E(^fe+i(i) - X k {t)) ■ P{k < n R {t)} 



k=0 



1 E(K(t)|)= r 



g — A g — A 

Similar considerations give the result for L(t). □ 
Lemma 3.6. 

Var(J ye > A (i)) < XE(\[V*t] - Q a (t)\) + 2t(p - q)X(g - A) + C{g - A)" 1 

Proof. The variance Var in the statement is taken in the five-process coupling 
where the Bernoulli distribution n is initially perturbed at the origin. Denote 
by Var A variance in the stationary process with initial invariant distribution /i^ . 
By the conditional variance formula 

Var A (J^ A (t)) = E A Var A (J ye < A (t) | r/ (0)) 

+ Var A E A (J^' A (<)M0)) 
>(l-A)-Var A (J l/e ' A (i)ho(0) = 0) 



= (1 — A) • Var( J 1/e,A (t)). 

Apply Proposition 12.11 to the first term Var A ( J yQ > x (t)). The conditional ex- 
pectations on the right-hand side of (|2.6[) match exactly the marginals of the 
five-process coupling constructed earlier in this Section, and so we come back to 
the present setting: 

V a r(J ve -\t)) < -±-Var x (J ve *(t)) = ^ A(l - \)E(\{VH] - Q(t)\) 
< XE(\[V e t] - Q a {t)\) + XE(\Q(t) - Q a (t)\). 
Write the last term above as 

(3.15) E(|Q(*) - Q a (t)\) = E(Q(t)) - E(Q B (t)) + 2E([Q(t) - Q a (t)]~). 
Proposition 12.11 gives 

E(Q(f)) - E(Q„(t)) = 2t(p - q)(g - A). 
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The last term in ([3~T5|) is treated separately for TASEP (p = 1) and ASEP 
(0 < q = 1 -p < p < 1). 

Consider first the totally asymmetric case p = 1. Initially Q(0) = Q a (0) = 0, 
and we show Q(t) > Q a {t) for all t > 0. The following situations cover all cases 
where Q and Q a could interchange positions: 

• If Q(t) = Q a (t) = i, then at site i (|3.5[) has the unique solution %(t) = 
CL^(t) = and r]f(t) = LOi{t) = 1. A right step of Q a without Q would 
require a brick on column i of lo, but not of r\ + . This is impossible by the 
basic coupling and T]f +1 {t) = rn+i(t) < uii + i(t). A left step of Q without 
Q a would require a brick on column i — 1 of r), but not of uf~ . This is again 
impossible by ?7,_i(t) < = cj i l 1 (t). 

• If Qa(t) = i and Q(t) = i + 1, then (|3.5|) has the unique solution ?fo(t) = 
r)f(t) = ur(t) = 0, Wi (t) - 1, Tft+iW = 0, and ^+ x (t) = ^^t) = 

= 1. None of the processes can have column i grow in this situation, 
hence Q a and Q cannot interchange positions. 

We conclude that the second term on the right of (|3 . 1 5|) is zero in the totally 
asymmetric case. 

For ASEP we use conditional coupling with event B of (|3.7[) . Let B = 
(PjS})" 1 < oo. 

E([Q(t) - Qa(t)}-) = E([Q(t) - Q a (t)}- | B) < E([L(t) - R(t)]- \ B) 

< BE([L(t) - R(t)}-) < BE(\R(t) - X(t)\) + BE(\X(t) - L(t)\) 

by the previous lemma. 

Collecting terms completes the proof of the lemma. □ 

Remark. The constant B = (P{Z?}) _1 diverges as p /* 1 and hence the same 
is true of C in the last term of the bound in Lemma \WM This suggests that there 
should be a more effective way to bound the variance of J v ' X (t) of the A-system 
in terms of quantities computed in the p-system. 

Lemma 3.7. 

Var(J e (i) - J ye > x (t)) < 2E(\[V e t] - Q a (t)\) + 4t(p - q)X(g - A) + C(g - A)" 1 . 
Proof. Similarly to the previous proof, we write 

Var e (J e (f)) > g- Var(J e (t)). 

Then we proceed by 

Var(J e (t) - J ye ' x (t)) < 2Var(J e (t)) + 2Var( J v ^ x (t)) 

< 2g- 1 Var e (J e (t)) + 2Var( J ye < X (t)) 

< 2(1 - g + A)E(|[^] - Q„(t)|) + 4t(p - q)X(g - A) + C(g - A)" 1 
utilizing the previous lemma. □ 
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We come to the lemma that summarizes all the previous estimation. 



Lemma 3.8. For any real u > 1 and time t > 0, 

(3.16) P{Q a (t) > Au+[VH]} < C 3 ^ ■ E(|Q„(t) - [VH}\) 

t 2 t 3 c ifl ~\ 

+ C 4 ^ + C 5 ^ + d expj -C 2 -j \ + e~ 2u . 
u 6 u° I t z ) 

Proof. Set b = 2(7 A g)(p — q) where 7 is the constant from Lemma 13.31 We 
proceed by cases. 

Case 1: 1 < u < bt. Suppose first u is an integer as it was in the proof of 
(13. 9[) . Throughout density g has been fixed, and now we also fix 

X = g 



2t(p - q) ■ 



The constraint on u guarantees that A > and g — A < 7 which was required 
for (|3.9p . The point of this choice of A is to maximize the lower bound inside 
the probability on the right-hand side of ()3.9|) . So, continuing from (|3.9|) with 
Chebyshev's inequality and the previous lemma, 

P{Q a {t)>Au+[VH]} 

< AP{*(t) - J^(t) > ^- } } + C X cx P {-c4} 

< C 3 ^Var(J^(i) - J VS '\t)) + C\ exp{-C 2 ^} 

t 2 t 2 t 3 ( II 3 1 

< C 3 — • E(|Q (t) - [V e t]\) + C A - + C 5 - + d exp{ -C2-3-}. 

f \Jb '\Xi 1-k ^* V * 



Extension from integral u to real u is achieved by adjusting constants on the last 
line above. 

Case 2: bt < u < 2t. Then bu/2 < bt and 

P{Q a (t) >Au + [VH]} < P{Q a (t) > 4 • bu/2 + [VH]}. 

Case 1 can be applied with u replaced by bu/2, at the price of adjusting some 
constants with powers of b. 

Case 3: u > 2t. Since Q a {t) is bounded above by a rate one Poisson process, 

P{Q a (t) >4u+ [VH]} < P{Q a (t) > 3u} < e~ 2u 

for all t. Combining the bounds from the cases proves the lemma. □ 



We are ready to complete the proof of the upper bound of Theorem 12.2 
Abbreviate 

^(t)=E(\Q a (t)-[VH}\). 
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First we need to complement p,16[) with a matching lower tail bound for 
Q a (t). This can be derived by arguments analogous to the ones we have pursued 
throughout Section O 

Alternately, we can derive the lower tail bound by a particle-hole interchange 
followed by reflection of the lattice. Define u3,(i) = 1— 0Ji(t), an ASEP with density 
g = 1 — q and rightward jump rate p = 1— p. Its second class particle has position 
Q(t) = Q a (t). Let (•) be the process obtained from £)(•) through a reflection 
about the origin. Then Q_ R (-) is ASEP with the original parameters (p, q). The 
second class particle of this process is at position Q(t) R — —Q{t) = —Q a {t), 
and the characteristic speed is V e,R — (p — q)(l — 2(1 — g)) = —V e . Since 
[-VH] = -~[VH], the expectation 

¥(t) - E(|Q a (i) - [VH]\) = E{\Q(t) R - [V^ R t}\) 

is the same for u>(-) and Q_ R (-)- Consequently 

P{Q a (t) - [VH] < -Au} = P{Q{t) R - [V s ' R t] > Au} 

which is again bounded by the right-hand side of (|3.16[) . 

Introduce a large constant 4 < r < oo. Consider t > 1 and u > rt 2 l z . Then 
we can combine the upper and lower tail bounds whose common right-hand side 
is given in (|3.16p . replace Au by u (we made sure u/A > 1), and simplify a little 
to arrive at 

P{\Q a (t)-[VH]\>u} 

( 3 - 17 ) i 2 / t 2 r « i \ 

< Cl - 4 nt) + C 2 (r) ( + ex P {-C 3 (r) Ws } ) . 

Constants were renamed and their dependence on r expressed in the notation 
because now it is of importance that C\ in front of ^S(t) does not depend on r. In 
particular, the exponentials in (|3.16| were combined via e~ 2u < exp(—2ut~ 2 ^ 3 ) 
and exp(-C 2 uH- 2 ) < exp(-C 2 r 2 ut- 2 / 3 ). 

Let 1 < m < 3. Integrate bound (|3.17|) over u S [rt 2 ^ 3 , oo): 

E(\Q a (t)-[VH]r) 

OO 

(3.18) < r m t 2m / 3 +m J P{ \Q a (t) - [VH}\ > uju^ 1 du 

rt 2/3 

< Cir m-4^ (t)t 2+(2/3)(m-4) + 

Constants were renamed again and their dependence on m ignored. To get the 
final bounds, take first m — 1 in (|3.18| to get 

#(t) < C ir - 3 y{t) + C A (r)t 2 ' 3 . 

Since C\ is independent of r, fixing r large enough gives ^(i) < C^{r)t 2 / 3 . Putting 
this bound back on the last line of (|3.18p then gives 

E(\Q a (t)-[VH]\ m ) <C 6 (r)t 2m ' 3 

for 1 < m < 3. The upper bound of Theorem 12.21 has been proved. 



20 



4 Lower bound 



The lower bound is proved by perturbing an initial equilibrium on a segment 
of the lattice. Again we begin with a description and some properties of the 
coupling. 

Perturbing a segment initially 

Recall again the characteristic speeds V e = (p—q)(l—2g) and V x = (p—q)(l—2X). 
We assume g > A, hence V e < V x . Throughout this section u > denotes a 
fixed positive integer, and 

n = [V x t] - [V e t] + u. 

To begin define an initial product distribution on two configurations (77(0), £(0)) 
by describing the marginals on each lattice site i: 

{ (77,(0), Ci(0)) ~ M of (03D if i < -n, 

r?i(0) = 0, 0(0) ~ Me ifi = -n, 

?7i(0) = Ci(0) ~ Ma if -n<i<0, 

IMO), Ci(0)) ~ m °f (ED if * > 0. 

Let Q(~ n \t) be the position at time f of a second class particle in the process 
rj{ ), initially at site Q^~ n \0) = —n. Note that site — n was set vacant for 77(0). 
As before rj(t) together with the particle at Q^~ n ^(t) make up the process r] + (t). 
Except for this perturbation 77(0) starts in the Bernoulli fi distribution. The 
process ((■) initially has distribution /x^, except at sites {—n +1, . . . , 0} where 
the parameter g has been replaced by A. 
Define a third initial configuration by 

6(0) = 

We let all these processes evolve jointly in the basic coupling. The following 
majorizations are true initially and are preserved by the evolution: 

V{t) < i(t) < C(f) an d h C (t) < ht(t) 

where the last inequality is for column heights. 

Let us denote the net particle currents (|2.5p by J v,v and J V '^ in the respective 
processes rj(-) and The first observation is that Q( -n ) gives one-sided control 
over the difference of these currents. 

Lemma 4.1. There is a nonnegative process N(t) with constant Geometric(qfp) 
time marginals P[N(t) = k] = (1 — q/p)(q/p) k and such that for any V 6 R 

Q [ ~ n) {t) < [Vt] implies J v ^{t) - J v ^{t) < N(t). 
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Proof. Denote the positions of the ^ — 77 second class particles at time t by 

• • • < Y k (t) <■■■< y_ 2 (i) < F_i(0) < Y (t). 

As argued before, we can arrange for these particles to keep their labels, and then 
the ordering is preserved. Initially Yo(0) < —n = Q( _n )(0). 

Let m,Q(t) — max{fc : Yu{t) < Q^~ n \t)} be the label of the £ — rj particle at or 
closest to the left of Q^ n )(t). Initially m Q (0) = 0. Once Q^ n \t) e {Y k (t)}, this 
containment property will hold forever because the basic coupling will preserve 
the ordering ry + < £ if this is ever established. 

We claim that ijiq remains zero while Q^~ n ^ is disjoint from the second class 
particles {F:}. This follows from showing that there is no jump which swaps the 
ordering F) < Q' - ™'- If Yo = i and Q^~ n ^ = i + 1, then we have = rji = r)f < 
£i = 1 and = r}i+i = Ci+i < Vi+i ~ 1- ^ n this situation column i of £ can 
increase, or the same column of r/ + can decrease. Neither of these steps can swap 
the positions of Fo and Q^ n \ instead they make Fo = Q^ n%> ■ 

Once Q(~ nS> is riding on the {Fc} particles so that actually 

Q(- n \t) = Y mQ(t) (t), 
its label rag evolves in the basic coupling as follows. 

• Suppose Q' - "- 1 = Fc = i and Fi+i = * + !•• A right Poisson arrow (i — + 
appears at rate p and increases ttiq from fc to fc + 1. 

• Suppose Q^~ nS> =Yk = i and Fc_i = i — 1. A left Poisson arrow (£ — ► i — 1) 
appears at rate g and decreases toq from fc to fc — 1. 

When f — 77 particle F mQ ( t ) itself jumps, Q^~ n > jumps with it. 

To get bounds on Q^ n \t) we introduce a suitable steady-state object. As we 
did with the £-marks (|3.3[) for the upper bound, we introduce a further classifica- 
tion among the £ — rj second class particles so that exactly one of them has priority 
over all the others. This special particle is marked by the label m(t). The Poisson 
arrows move m(t) exactly the same way as rriQ(t) on the labels {—00 < fc < 0}: 

• If Y m m(t) — i and F m ( t ) +1 (t) = £+ 1 and there is a right Poisson arrow 
(i — > i + 1) then m(£) increases by 1. This happens at rate p. 

• If Y m r t )_x(t) = i — 1 and F m ( t )(i) = i and there is a left Poisson arrow 
(i — > i — 1) then m(t) decreases by 1. This happens at rate g. 

From these rates we see that m(t) behaves like a birth and death chain on Z_ 
whenever adjacency of £ — 77 second class particles permits m(t) to jump. Without 
obstruction this birth and death chain would have reversible measure n(k) — 
(1 — q/p)(q/p)^ for fc < 0. We give m(0) initial distribution n. Then the 
argument given for Proposition 13.11 implies that for each fixed time t we have 
P[m(t) = fc] = 7r(fc). 

We have arranged m(0) < mg(0) initially. Since Q(~ n ^ cannot hop over Yq 
without joining it, the identical responses of m and vhq to the Poisson arrows 
implies that m(t) < mg(i) for all time. 
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To connect with currents, note that -mg(t) equals the number of £ — 7] 
second class particles strictly to the right of Q(~ n \t) at time t. Since these 
particles started off in (—00, 0] and they account for all the discrepancies between 
the processes £ and 77, this number equals the height difference h^_ n) ^(t) — 



hi 



Q(-")(t) 



(t). Thus under {Q ( " n) (» < [Vt]} we have 
-m(t)> - m Q {t) = h s Q( _ n)(t) (t) - h v Q( _ n){t) (t) 



hi vt M-h*{t) + 



i=Q(-«)(t) + l 



E 



>ftfvt](*)- ft Fvt](*) 

>^t](i)-^ t] (i) = ^' C (i)-^)- 
To obtain the statement of the lemma take N(t) = —m(t). 



□ 



Let w(-) be a stationary ASEP started from the [i Bernoulli distribution. 
The next lemma gives a way to compare the distributions of Q and w. 

Lemma 4.2. Denote by P w and P'' the probability of events that depend only on 
the respective processes w(-) and ((■). Then 



P C (-) < P w (-)^ -exp 



n(g - A) 5 



2e(i - q) 



Proof. Let 



^= E 6(0). 

i— — 7i+l 

Z has a Binomial(n, A) distribution ^ A (z) = (")A Z (1 — A)"~ z for 2 = 0, . . . , n. 
We use the Cauchy-Schwarz inequality below to perform a change of measure on 
this binomial distribution. The binomial mass functions in the second line are 
easily added up. 



P<(.) = ]TP c (-|Z = z) [v«(z)]i 



z=0 



< 



.[E[ pC (-^ 

n 

= [E[ pC (-^ 

n 

<[E pC (-^ 

z=0 



= Z 



= Z 



)] 2 i/«(*)] 
)] 2 */>(z)] 

z)u Q {z) 



v x {z) 
[ve(z)]i 

1 n r \ 



(*)] 



2 1 



1 + 



(g-X) 



2 , 



cxp 



g(l - g) 
n(g - A) 2n 



2e(i - e) 
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All that is left is to recognize that P^( • | Z = z) is the probability of a process £(•) 
whose initial distribution is Bernoulli^) outside {— n + 1 . . .0}, with z particles 
distributed in that interval with each configuration equally likely. Summing these 
with the Binomial(n, g) coefficients v e (z) gives the Bernoulli initial distribu- 
tion of the process w(-). □ 



Proof of the lower bound 

With these preparations we are ready to prove the lower bound. As for the upper 
bound, Q a (t) is the position of a second class antiparticle started from the origin 
on a ^-equilibrium process w(-) initially perturbed by setting u>o{0) = 1. Our 

target quantity is abbreviated as before by 'i'(i) = E(\Q a (t) — [V^ e i]Q. We start 
with bounding the probability of the complement of the event in Lemma |4. II As 
before u is an arbitrary but fixed positive integer and n — [V t] — [V e t] + u. 

Lemma 4.3. 

P W H0(t) > [VH]} < m + afr-gXe-A)-^ + C 

u u u{g — X) 

Proof. For this proof, Q(t) will refer to a second class particle started from the 
origin on a process fy(-) in fi^ distribution, except for 770(0) = 0. Translation 
invariance implies 

P{Q { ~ n) (t) > [V g t]} = P{Q(t) - [V x t] >u}< E(IQW ~ [yAt]l) 

< E(|Q(t) - Q a (t)\) | E(|Q a (t) - [VH]\) | [VH] - [VH] 
u u u 

where we introduced the couplings of Section [3] As in Lemma l3.6[ the first term 
is bounded from above by 2t(p — q)(g — \)/u + C/[u(g — A)]. The second term is 
^f(t)/u, and the third term is bounded from above by 2t{p~q)(g — X)/u + 2/u. □ 

Lemma 4.4. For any < K < (p — q)t(g — A) 2 , 



(p-q)t(g-\y-K ^l2g(l-g). 
4A*(t) 8t(p — q)X(g - A) C C 



K 2 K 2 K 2 (g-X) K-8' 

Proof. Lemma SH] leads to 

P{Q(-")(i) < [VH]} < P{J ve <t(t) - J ve ^{t) < N(t)} 
(4.1) < P{ J ye ^ (t)<K+(p- q)t{2gX - A 2 ) - 1} 

K 



(4.2) + P{ J ve '\t) - (p - q)t(2gX - A 2 ) > - + 3} 

(4.3) +P | iV (t)>^-4}. 
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We apply Lemma 



lity: 



to line (|4.1[) to bound it by the equilibrium /i -probabi- 



[P e { J e {t) <K+{p- q)t(2g\ - A 2 ) - 1}] 2 • exp 

< [P e { J e (t) <K-(p- q)t(g - A) 2 }] * • exp 
[Var(je(t))]V2 \n{g-\) 2 



n{ Q -Xf 



2q{1 - q) 
n{g-\) 2 



2g(l - g) 



< 



(p - q)t(g - A) 2 - K 



(p-q)t(g-\y 



K 



■ exp 
• exp 



2g(l-g). 
n{g~\) 2 



2f?(l - g) 



The term —1 was subsumed in errors caused by integer parts when we centered 
J«(t). 

A simple coupling consideration shows that E(J y ' ,T, {t)) differs by at most 
one from the same expectation taken under an unperturbed fj, initial condition. 
Thus taking integer parts again into account, line (|4.2|) is bounded by 



P{j^(i) > \ } 



K\ 4Var(J 



v e . 



K 2 



Lemma 13.61 can be applied to bound this variance even though the second class 
particle now starts at — n rather than at the origin. The only change needed in 
the proof of Lemma 13.61 is in the calculation (I3.14|) where one must condition on 
i]- n (0) instead of on r?o(0). Hence we can continue from above to bound line (14. 2p 
with 

4Att(t) 8t(p-g)X(g-X) C 
K 2 K 2 K 2 {g-\y 

Lastly, the geometric probability (|4.3|) is bounded by CJ (K ■ 
shev's inequality. 



by Cheby- 
□ 



Now the last step of the lower bound of Theorem 12. 21 By Holder's or Jensen's 
inequality it suffices to prove the case m — 1, in other words that 

liminf t~ 2/3 ^(t) > 0. 

t — >oo 

In the two last lemmas take 

u = \ht 2/3 ] , g-X = bt- 1/3 , and K = bt 1/3 , 

where h and b are large, in particular b large enough to have b < (p — q)b 2 so that 
K satisfies Lemma FOl Then 



[V x t] 



V e t] + u< [2(p - q)b + h]t 2 ' 3 + 3 < Ct 2/3 

and Lemma 



for large enough t. We can simplify the outcomes of Lemma 
to the inequalities 



(4.4) 



P{Q(-")(i) > [VH]} < C 



£2/3 



Ab 
~h 



C 
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and 



(4.5) PW C-n) (t) < [V H ]} < C(^lf 2 + C§± + 1 + ^. 

The new constant C depends on b and h. 

The lower bound now follows because the left-hand sides of (|4.4p ~ (|4.5p add up 
to one for each fixed t, while we can fix b large enough and then h large enough 
so that Ab/h + 8/6 < 1. Then t~ 2 / z ^(t) must have a positive lower bound for all 
large enough t. This completes the proof of Theorem 12.21 

Remark. The key properties, namely /" ^ for the flux and asymmetry of 
the jump kernel, were used somewhat surreptitiously in the above proof. As in 
calculation (|3.13|) . the nonvanishing of /" renders E[J e (t) — J ye ' r, {t)] of order 
i 1 / 3 . This allowed us to take K of order t 1 ^ 3 , then g — A of order i -1 / 3 , and u of 
order t 2 ^ 3 . This way the final bounds (|4.4[) - (|4.5[) have ^(t) divided by the correct 
order t 2 / 3 . 

Asymmetry also came in: if p = q = 1/2 then the means of the equilibrium 
currents are zero and there would be no deviations to take advantage of in the 
proof of Lemma (|4.4[) . 
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